Interaction between superconductive films and magnetic 

nanostructures 



Lars Egil Helseth 

Department of Physics, University of Oslo, P.O. Box 104-8 Blindern, N-0316 Oslo, Norway 



Abstract 

The interplay between magnetic and superconductive films is studied. The 
generalized London equation for this system is solved, and the magnetic fields, 
the energy and the interaction forces are computed. In particular, we focus 
on how to manipulate vortices using magnetic nanostructures. 
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I. INTRODUCTION 



The structure of vortices in thin films was first investigated in detail by Pearl [l]. He 
found that such vortices interacts mainly via their stray field, which extends far into the 
nonsuperconducting medium. Later, this theory has been extended to thin film systems 
with and without anisotropy It has also been of cnsiderable interest to understand 

Josephson vortices in thin film junctions. To this end, a very interesting result was obtained 
by Kogan et al. who found that the field from a Josephson junction is a superposition 
of fields from Pearl vortices along the junction with a certain line density. This result was 
derived by noting that the Pearl solution is the Green's function for arbitrary sources. 

The possibility of enhancing the critical currents by magnetic nanostructures have been 
studied by many researchers ||[7|]. In particular, Bulaevskii et al. showed that magnetic 
domain structures in a magnetic film in close contact with a superconducting film may en- 
hance pinning of vortices, since this gives an opportunity to pin the magnetic flux of the 
vortex rather than its core J7J. It was suggested that the pinning of vortices in supercon- 
ductor/ferromagnetic multilayers can be many times greater than the pinning by columnar 
defects. 

In general, the interaction between superconductivity and magnetism has been studied 
for several decades. Systems composed of magnetic and superconductive materials are of 
interest not only because they are model systems for the interplay of competing supercon- 
ducting and magnetic order parameters, but also because of numerous possible applications. 
Recently, the development of magnetic thin film technology has triggered a lot of interest 
in this field. Of particular importance has been the possibility of examining the interaction 
between superconductivity and magnetism in high-temperature superconductors P-TT5[|. 

In a recent paper we estimated the interaction between a bulk superconductor and a Bloch 
wall fl6|]. It was shown that the Bloch wall shrinks on cooling through the superconducting 
transition temperature. Furthermore, we suggested that if the magnetization of the Bloch 
wall is strong enough, it may capture and move the vortices. 

The aim of the current paper is to study the interaction between Pearl vortices and 
magnetic nanostructures in more detail. We solve the generalized London equation, and 
compute the magnetic fields, the energy and the interaction forces. In the case of a movable 
magnetic nanostructure (e.g. a domain wall), one may be able to manipulate the static as 
well as the dynamic behaviour of the vortices. Thus, tunable and movable nanomagnets 
may be used as vortex manipulators. 



II. THE THIN FILM SYSTEM 

Consider a thin superconductive film and thin magnetic film, both located at z=0 with 
thicknesses much smaller than the penetration depth of the superconductor, see Fig. [I]. The 
two are separated by a very thin oxide layer (of thickness t) to avoid spin diffusion and 
proximity effects (which may supress the superconductive and magnetic order parameters) 
|T7| |. In general, the current density is a sum of the supercurrents and magnetically induced 



currents, which can be expressed through the generalized London equation as 



V x J = —H + ^V(p)e z + V x V x M v , (1) 
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where J is the current density, A is the penetration depth, H is the magnetic field and 
V x M v is the magnetically induced current. The vortex is aligned in the z direction, 
and its source function V(p) is assumed to be rotational symmetric. In the case of a Pearl 
vortex we may set V(p) = (§ / fi )5(p), where $ is the A ux quantum and /i the permeability 
of vacuum. The requirement of rotational symmetry should be relaxed when dealing with 
certain types of anisotropic superconductors, which are known to have azimuthal corrections 
to the vortex core [ IE \ . In this case one should also use the anisotropic version of Eq. ([I]) . 



Recently, the Ginzburg-Landau theory was used to calculate the vortex states in a thin 
mesoscopic disk in the presence of step-like external magnetic field [J19||. Naturally, one 



may expect more accurate results from the Ginzburg-Landau theory, in particular when the 
dimensions are small. On the other hand, the generalized London theory considered here 
gives significantly reduction in computational efforts, and some important points could be 
deduced without numerical analysis. Moreover, the inherent shortcomings of the London 
model can be relaxed by introducing corrections to the vortex core. Thus, we believe that 
the generalized London model adopted in this paper is a reasonable first approach. 

We solve the current London equation by adopting a generalization of the very useful 
approach of Refs. |5|,|20[|. The current can only flow in a thin layer of thickness d (d < A), 



and it is therefore necessary to average over this thickness and consider the z component 
only 

/ d I s f)T s \ 

H z + X e h^-TT 1 )=V(p) + X e (V XV X M) z , (2) 



dx dy 

where A e = \ 2 /d is the effective penetration depth, and M = dM v . The sheet current 
flowing in the thin layer is now given by J s = dJ. The Maxwell equation V X H = J gives 

j = d_H^ _ (My j = dH* _ dlh (3) 

x dy dz v dz dx 

Since all derivatives d/dz are large compared to the tangential d/dp, we may set 

J X ~ H y — H y , jy ~ " , 

where Hf and (i=x, y) are the components at the upper and lower surfaces, respectively. 
Since the enviroments of the upper and lower half-spaces are identical, we have = —Hf, 
which results in 

J s x « -2H+, J s y w 2H+ . (5) 

Here we will only consider the upper half-space. Using V • H = 0, Eq. (H) becomes 

dFf 

H z - 2X e -^f = V(p) + A e (V x V x M) z . (6) 

We now apply the superposition principle to separate the contributions from the vortex and 
the magnetic nanostructure. The vortex part of Eq. (||) can be written as 

d FT 

H vz - 2X e ^ = V(p) , (7) 

whereas the magnetic part is 

dM 

H mz - 2A e ^^ = A e (V x V x M) z . (8) 
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A. Vortex solution 



In order to solve Eq. ([/]), it is useful to note that Vxi3" = V- i3" = outside the thin 
films. Therefore, we can introduce a scalar potential <ft which vanishes at z —>■ oo 

"I poo poo 

MP: z ) = /_ oo y_ oo ex p( ifc • p - fcMKfc , (9) 

where k = ^jk\ + and i/ = V0. It is helpful to note that H z {k) = —k(j){k) for the 
upper half-space. 

Applying the 2D Fourier transform to Eq. ([7]) we obtain 

« fc > = ' (10) 

In the case that the vortex is displaced a distance p from the origin we may set 

, v V(k) exp(-ik • p ) 

= (l + 2\ e k)k • (H) 

However, for the moment we assume that the vortex is located at the origin, and therefore 
the resulting scalar potential is 



OO POO 



Due to the rotational symmetry the potential is found to be 

*- (ft '» = -sr v(i 'i^" 5 (- w *' (i3 » 

where we have used that 

/ exp(ikpcos<fi)d(J) = 2irJ (kp) . (14) 
Jo 

To obtain the magnetic field components (in the radial and z direction), we apply the 
following formula 

^-Ukp) = -kJ^kp) , (15) 
dp 



and find 



H vz (p,z) = ±- [°° kV(k)-^^-ex V (-k\z\)dk , (16) 

h '^ z) = i r kv ^T^ e ^- k ^ dk • < i7) 

In the case of a Pearl vortex the magnetic potential and fields reduce to 
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B. Magnetic solution 

In general, Eq. (H) can be solved for arbitrary magnetization distributions. For simplicity, 
we will here consider only magnetization distributions which are directed in the z direction 
and have no volume charges, V • M = 0. Inserted in Eq. (||) this gives 

OH 

H z - 2A e ^ = -A e (V 2 M) 2 . (21) 

The potential from the magnetic nanostructure can be obtained in the same way as for 
the vortices. We will here only work out explicitly the solution for two-dimensional (2D) 
magnetic distributions, but the ID solution is found in a similar manner. Applying the 2D 
Fourier transform to Eq. ([5]) results in 

« fc > = • (22) 

which gives the following magnetic potential: 

A e [°° f°° kM{k) ^ 2 



m (p, z) = -7^31 y_ oo y_ oo 1 + 2 \ k ex v( ik ' P) ex P(~ k \z\)d k . (23) 
In the case of a rotational symmetric magnetization distribution, we may write 

m (p,z) = r J ( k p)^&-eM-k\z\)dk . (24) 
Zn Jo 1 + Z\ e k 

Using that H = V0, we can calculate the radial and vertical field components: 

H mz (p,z) = ^ r J (kp)^^-ex V (-k\z\)dk , (25) 

Zll Jo 1 + Z\ e K 

and 

H mp (p, *) = ^ jH eM-k\z\)dk . (26) 

So far we have neglected the thickness of the oxide layer. However, to a first approximation, 
it is easy to show that the only effect of this layer is to shift the z coordinate, z — > z — t. 
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C. Energies and forces 



We now work out the energies and forces in the case of a circular symmetric 2D magne- 
tization distribution, but will later also discuss briefly a ID distribution as well. Consider a 
vortex located a distance po from the origin, and a magnetization distribution centered at 
the origin. The three contributions to the energy are the magnetic self-energy, the direct 
interaction energy and the energy associated with the vortex. 

The magnetic self-energy is given by 

I /"OO rOO 

E m = --p / M • H m d 2 p . (27) 

Ji J —oo J —oo 

This expression can be transformed by using the following formula: 

/oo/'OO 1 roo roo 

/ A(p)B(p)d 2 p = — / / A(—k)B(k)d 2 k . (28) 
-oo J —oo I Z7T J J —oo J —oo 

Thus, in the case of a circular symmetric magnetization distribution we have 

MP MWMW 
4tt Jo l + 2X e k ' V 7 

where we have used that M(-k)=M(k) and also that d 2 k —>■ kdkd<p (and integrated over the 
azimuthal angle). 

The direct interaction energy between the vortex and the magnetization is given by 

/oo roo 
/ M ■ H v d 2 p , (30) 
-oo J —oo 

where the vortex is assumed to be displaced a distance p from the origin. Then, in the 
cylindrical symmetric situation one finds 

The associated force can be found by taking the derivative with respect to po, resulting in 

F m = -!« r k 2 V(k) M{k)J f^ dk . (32) 
The energy associated with vortex field is expressed as 

E 3 = J v (^PoH 2 + \n s mv 2 )j dV , (33) 

where n s is the density of superconducting electrons and v s is their velocity. This expression 
is most easily transformed by applying the formalism of Ref . |2(J . In particular, the energy 
associated with a Pearl vortex can be evaluated to logarithmic accuracy as 

- %ln , (34, 



Aicp X e \ f 

where ^ is the coherence length. More accurate expressions are obtained by introducing the 
appropriate vortex source function, but this will not be considered here. 
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III. INTERACTION WITH A BUBBLE 



Consider first a magnetic bubble that can be moved, generated or annihilated by external 



stress patterns or magnetic fields pl| . Let the vortex be separated a distance po from 
the magnetic bubble located at the origin. The system has cylindrical symmetry, and the 
magnetic bubble have a magnetization vector pointing in the z direction given by 



M = 

The Fourier transform of this distribution is 



M e z if p < W 
if p > W 



M(k) = 2 7 rAW Jl( f y) , (35) 



which generates the following magnetic field components: 

k 2 J x {kW) 



H b z (p,z) = X e WM I" Joikp) * l Q , I exp(-k\z\)dk , (36) 
Jo 1 + l\ P k 



H b p (p,z) = \ e WM r Ji{kp) k l J ^} eM-k\z\)dk . (37) 
1 Jo 1 + 2A e k 



These components were also found in Refs. |I3|J14| using a different method. To obtain the 
total field one should add the vortex contribution. The energy is found to be 

$1 , /47rA e \ t ,^ rir f°° Ji(kW)J (kpo) „ , rir2 f°° , JlikW) „ 

(38) 

and the direct force between the bubble and the vortex is given by 

y^-^n mi ft. (39) 

Jo 1 + 2X e k 

Figure shows the force when W = A e , and we see that it falls off rather quickly. Moreover, 
we note that there are peaks at the maximum magnetization gradients, as expected. A 
vortex in the vicinity of the bubble will be either repelled or attracted to the center of the 
bubble, where it eventually comes to rest. The vortex-bubble system may still attract more 
vortices, but now vortex- vortex repulsion must be taken into account, which is outside the 
scope of the current study. 

Consider now an annular structure as illustrated in Fig. [| where the inner dot has 
magnetization Mi and the outer annulus M2. Then the force becomes 



-MiWi$ 



As and example, Fig. (| displays the force when Mi = IM2I, W\ = 0.5A e , W3 = 2Wi, 
and W2 = 1.5Wi. The solid line shows the case of a positive M2, whereas the dashed line 
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corresponds to negative M2. A slowly moving vortex approaching the annular structure 
from the outside will experience attraction if M 2 > 0, and its equilibrium position is directly 
under the annulus. If M 2 suddenly switches sign, then the vortex is attracted to the center 
of the annular structure (i.e. directly under the dot) and kept there. Thus, this annular 
structure could be used to catch and trap vortices by switching the sign of the outer annulus. 
This could in principle be done by using a soft magnetic material with at the outer annulus, 
whereas the inner dot consists of a hard magnetic material. Naturally, the magnetic forces 
must overcome the Lorentz forces generated by the external magnetic field. Moreover, it 
is probably not so easy to fabricate such a structure. Therefore, it is doubtful whether 
this circuit could find any practical realization, although it is a nice example of a vortex 
manipulation structure. 



IV. INTERACTION WITH A ID MAGNETIZATION 

Most domains and walls are not cylindrical symmetric. In fact, they are more often one 
dimensional (ID), and therefore the corresponding London equation should be solved for 
this symmetry. Consider a ID magnetization distribution with the vector pointing in the z 
direction. Then we may apply the methods presented in this article (using the ID Fourier 
transform), and find the following expression for the scalar potential: 

<l> m (x,z) = -— [ k * M ( k *)™ s ( k * x ) e x p (-k x \z\)dk x , (41) 
n Jo 1 + ZA e K x 

where we have assumed that the Fourier transform of the magnetization is real. The mag- 
netic field components are given by 

A e [<*> 2 M{k x )cos{k x x) 

H mz (x,z) = — k x — - • ~~~ ~ exp(-k x \z\)dk x , (42) 

n Jo 1 + Z\ e k x 

u ( \ Ae t°° U 2 M ( k x)sin(k x x) 

H mx (x,z) = — k x — — — — exp(-k x \z\)dk x . (43) 

it Jo 1 + 2X e k x 

The energy of the system is found to be 

E = _*°-l„ ( ^) - *± f M ^)™< k ^) dkx _ ^0 r kl M\k x ) ^ 
47r/i A e \ £ / 11 ^° 1 + 2A e &z In Jo 1 + 2A e fc x 

To find the interaction force between the magnetization distribution and the vortex, we 
simply take the negative derivative of the interaction energy, and obtain 

fm = -*h rk, M ^if;^ d K. (45) 

71 Jo 1 + 2X e k x 

As an example, Fig. |5| shows the z component of the magnetic field generated by the 
following magnetization distribution (which is looks like a Bloch wall separated by in-plane 
magnetized domains): 



S 



f M 5(z)e z if -W<x<W 
{ if |x| > W 

Here we have chosen W=A e /50, z=0 and t=A e /2000. Note in particular the peaks near 
the edges of the Bloch wall. In general, it is also seen that there is a small field directly 
above the magnetization distribution which is not screened by the superconductor. It can 
be shown that the corresponding force on the vortex has a shape similar to that of Fig. |2|, 
with its maximum near the edges of the magnetization distribution. 



V. CONCLUSION 

The interplay between magnetic and superconductive films is studied. The generalized 
London equation for this system is solved, and the magnetic fields, the energy and the 
interaction forces are computed. In particular, we focus on how to manipulate vortices using 
magnetic nanostructures. The formalism presented here could be generalized to arbitrary 
magnetization vectors, and we therefore hope that it will be useful in future studies. 
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FIGURES 

FIG. 1. The system under study. The thin oxide layer separating the the superconductive and 
magnetic films is not shown. 

FIG. 2. The force between a vortex and a bubble when W=A e . 

FIG. 3. An annular structure for manipulating of vortices. 

FIG. 4. The force on a vortex near the annular structure of Fig. |3|. 

FIG. 5. The magnetic field generated by a Bloch wall with width 0.02A e . Here we have assumed 
that z=0.0005A e . 
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Figure [| 
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Figure [| 
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